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Abstract 

A new minimal coupling method is introduced. A general dissipa- 
, tive quantum system is investigated consistently and systematically. 

Some coupling functions describing the interaction between the sys- 
tem and the environment are introduced. Based on coupling func- 
tions, some susceptibility functions are attributed to the environment 
explecitly. Transition probabilities relating the way energy flows from 
the system to the environment are calculated and the energy conserva- 
tion is explecitly examined. This new formalism is generalized to the 
dissipative scalar and vector field theories along the ideas developed 
for the quantum dissipative systems. 

Keywords: Dissipative Quantum Systems, Field Quantiza- 
tion, Environment, Scalar and Vector Fields, Coupling Func- 
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1 Introduction 

There are basically two approaches to study dissipative quantum systems. 
One is found in the interactions between two systems via an irreversible en- 
ergy flow PU |2] . The second approach is a phenomenological treatment under 
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the assumption of nonconservative forces |H1 0] . In studying nonconservative 
systems, it is essential to introduce a time dependent Hamiltonian which de- 
scribes the damped motion. Such a phenomenological approach for the study 
of dissipative quantum systems, especially a damped harmonic oscillator, has 
a rather long history. Caldirola and Kanai E] adopted the Hamiltonian 

H{t) = e-*£ + e*lmJZf, (1) 

which leads exactly to the classical equation of motion of a damped harmonic 
oscillator 

q + (3q + u 2 q = 0. (2) 

The quantum aspect of this model has been studied in a great amount of 
literature. In those studies some peculiarities of this model and some fea- 
tures of it have appeared to be ambiguous [Zj-jTH]- There are significant 
difficulties in obtaining the quantum mechanical solutions for the Caldirola 
-Kanai Hamiltonian. Quantization with this Hamiltonian violates the uncer- 
tainty relations, i.e., the uncertainty relations vanish as time goes to infinity 

Based on Caldirola-kanai Hamiltonian, it has been constructed equivalent 
theories by performing a quantum canonical transformation and has been 
used the path integral techniques to calculate the exact propagators of such 
theories, also the time evolution of given initial wave functions, have been 
studied using the obtained propagators [2*0] . 

In the framework of the phenomenological approach Lopez and Gonza- 
lez [21J have taken the external non conservative forces that has linear and 
quadratic dependence with respect to velocity. They have deduced classical 
constants of motion and Hamiltonian for these systems and eigenvalues of 
these constants have been quantized through perturbation theory. 

A simple pseudo-Hamiltonian formulation has been proposed by Kupriyanov 
and et al [22]- Starting from this pseudo- Hamiltonian formulation a con- 
sistent deformation quantization has been developed that involve a non- 
stationary star product and an extended operator of time derivative dif- 
ferentiating the star product. A complete consistent quantum- mechanical 
description for any linear dynamical system with or without dissipation has 
been constructed in this scheme. 

In another approach to quantum dissipative systems, one tries to bring about 
the dissipation as a results of an averaging over all the coordinates of the bath 
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system, where one considers the whole system as composed of two parts, our 
main system and the bath system which interacts with the main system and 
causes the dissipation of energy on it [2lI]-[3T)|. 

The macroscopic description of a quantum particle with passive dissipation 
and moving in an external potential v(x) is formulated in terms of Langevin- 
Schrodinger equation [2U E21 E3] 

mx + f dt'fx(t - t')'x{t') = -Vv(x) + F N (t). (3) 
Jo 

The coupling with the heat-bath in microscopic levels corresponds to two 
terms in macroscopic description. A mean force characterized by a memory 
function /j,(t) and an operator valued random force F/v(£). These two terms 
have a fluctuation-dissipation relation and both are required for a consistent 
quantum mechanical description of the particle. In [33], there are some 
models for interaction of the main system with the heat-bath which lead 
to the macroscopic Langevin-Schrodinger equation. 

The layout of the paper is as follows: In section 2, we consider a parti- 
cle moving in a one-dimensional external potential and take the absorptive 
environment of the particle as a Klein-Gordon field which interact with the 
momentum of the particle through a minimal coupling term. In this ap- 
proach, the Langevin-Schrodinger equation© is obtained as the macroscopic 
equation of motion of the particle and the noise force F N (t), is derived in 
terms of the coupling function and the ladder operators of the environment. 
By choosing a special form for the coupling function, a friction force propor- 
tional to the velocity of the particle is obtained . For an initially exited one 
dimensional damped harmonic oscillator, it is shown that the entire energy of 
the oscillator will be absorbed by its environment. In section 3, we are con- 
cerned with a particle moving in three-dimensional absorptive environment. 
In this section we model the environment by two independent Klein-Gordon 
fields B and B. The field B interacts with the momentum of the particle 
through a minimal coupling term and the field B, interacts with the position 
operator of the particle similar to a dipole interaction term. A generalized 
Langevin-Schrodinger equation is obtained which contains two memory func- 
tions describing the absorption of energy of the particle by its environment. 
As an simple example, the spontaneous decay constant and the shift fre- 
quency of a two- level system, embedded in an absorptive environment, is 
calculated. 
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In sections 4 and 5, the ideas set up in sections 2 and 3, are applied to 
quantum dissipative systems described by scalar or vector fields. 



2 One-dimensional quantum dissipative sys- 
tems 

2.1 Quantum dynamics 

Quantum mechanics of a one-dimensional damped system can be investi- 
gated in a systematic and consistent way by modeling the environment of 
the system with a quantum field B which interacts with the system through 
a minimal coupling term. For this purpose, let the damped system be a 
particle with mass m, under the influence of an external potential v(q). We 
take the total Hamiltonian of the system and it's environment as (3T] 

H=b^ + v{q)+H B , (4) 

where q and p are position and canonical conjugate momentum operators of 
the particle respectively and satisfy the canonical commutation rule 

[q,p]=ih. (5) 

In Hamiltonian (j3J), Hb is the Hamiltonian of the environment. Now we 
model the environment by a massless Klein-Gordon field, such that we can 
write the environment Hamiltonian as 

/+oo 
d 3 fcftu,,4(%(i), (6) 
-oo 

where cug, is the dispersion relation of the environment. The annihilation 
and creation operators 6g, oi, in any instant of time, satisfy the following 
commutation relations 

h(t),bl,(t)] = 6(k-k'). (7) 

Operator R in relation (J3J), has the basic role in interaction between the 
system and it's environment and is defined by 

R(t) = / d'kifi^Mt) + r(u % )bUt)), (8) 
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where the function f(oo^) is called the coupling function between the parti- 
cle and its environment. In definition of operator R, the following physical 
assumptions are considered: 

1. Since the total Hamiltonian (jlj is Hermitian, so the operator R should be 
a Hermitian operator. 

2. We are interested in linear dissipative systems and for such systems, 
the operator R should be a linear combination of annihilation and creation 
operators bj: and £>t. 

3. Operator R, should depend on macroscopic characters of the environ- 
ment, because we are interested in macroscopic equations of motion of the 
particle. It is remarkable to note that the macroscopic characters of the en- 
vironment are reflected in its macroscopic susceptibility against the motion 
of the particle. In the following, we see that the susceptibility, is related to 
both the dispersion relation and the coupling function f(u>t) and these 
considerations are contained in definition of the operator R. 

The Heisenberg equations for the position and momentum of the particle 

are 



where H is the total Hamiltonian (fl|). By eliminating p from relations (JHJ), 
we obtain 



i 



p- R 



rn 

dv 



(9) 



dq 




(10) 



Also using ((7j), the Heisenberg equation for b%, is 




i 



i 



(11) 



which formally can be solved as 



Hit) = b n (0)e-^ + jf dt'e-^ 



(12) 
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Now substituting b^(t) from this equation into (jSJ), we obtain 

R(t) = R N (t) + [ ]t] dt' X (\t\ - t')q(±t'), (13) 
Jo 

where the upper(lower) sign corresponds to t > 0(t < 0), respectively. In- 
spired by electrodynamics in a media, let us call the memory function 

S77- roo _, 

X (t) = - d\k\\k\ 2 \ fiu^ sin u,:t t >0 



h 

X (t) = t < (14) 

as the susceptibility of the environment. The operator Rn 

/+00 
^ d 3 k[f(^)b u (0)e-^ + r(^)4(0)e^], (15) 

is a noise operator. 

The noise forces are necessary for dealing with dissipative quantum sys- 
tems in a consistent way, without theses noises, some fundamental postulates 
of quantum mechanics like canonical commutation relations or Heisenberg 
uncertainty relations usually are violated. 

The form of the equation (fT3|) suggests that take this equation as a consti- 
tutive equation of the environment and interpret the operator R as some kind 
of the polarization of the environment, induced by the particle. A feature of 
this approach is its flexibility to choosing an appropriate dispersion relation 
ujj: and finding the corresponding coupling function for a given susceptibility. 
In fact, in definition (|14|) of the susceptibility, there are two main functions, 
namely the coupling function f(u)k) an d the dispersion relation u^, if one 
takes the simplest assumption uo^ = c\k\ for the dispersion relation and write 



as 



X (t) = 0, t<0, (16) 

then by inverting (|16|h the corresponding coupling function for a given sus- 
ceptibility, can be obtained as 

fl(y^ roo 

l/HI 2 = -^ji J o dt X {t) smut, u>0, 

|/H| 2 = 0, u = 0. (17) 
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Other choices for the dispersion relation, just lead to a redefinition of the 
coupling function and also more difficult mathematical expressions which 
basically may not allow to take an inverse of (fPfj) and find the corresponding 
coupling function in terms of the susceptibility x(£). It is remarkable to note 
that one can take the real dispersion relation of the medium and enrich the 
model in this way but in this approach this is only a choice. From now on, 
we take a linear dispersion relation as Ur — c\k\. 

Substituting R from (fT3]) into (fTUj) . one can show easily that the equation 
of motion for q is 

mq±jJ o dt' X (\t\-t')q(±t') = -—-R N , (18) 

this is the macroscopic Langevin-Schrodinger equation for a dissipative sys- 
tem [33] . If we take the Langevin-Schrodinger equation as the basic macro- 
scopic equation of quantum dissipative systems, then the microscopic inter- 
action of the particle with its environment can be macroscopically formulated 
by modeling the environment with a quantum field B and using the Hamilto- 
nian (j3J), which finally leads to the correct macroscopic Langevin-Schrodinger 
equation. 

In fact, the present approach can be considered as a generalization of the 
Caldeira-Legget model [23 , where for dissipative quantum systems, they 
model the environment of the system by a collection of harmonic oscillators 
such that the position operator of the main system, is coupled with all of 
the position operators of the environment oscillators. While the properties 
of the environment may in some cases be chosen on the basis of a micro- 
scopic model, this does not have to be the Caldeira-Legget model. As an 
example we mention an Ohmic resistor which as a linear electric element 
should be well described by the Caldeira-Leggett model. On the other hand 
the underlying mechanism leading to dissipation in a resistor may be more 
complicated than that implied by the model of a collection of harmonic oscil- 
lators. Modeling the environment with a collection of Harmonic oscillators, 
with a suitable interaction, means that an averaging over complicated mi- 
croscopic interactions of the particle with its environment is done and the 
environment is effectively equivalent to this modeling. 
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2.2 Quantum damped harmonic oscillator 

For a one dimensional harmonic oscillator with mass m and frequency a>o, 
we have v(q) = \muj 2 q 2 and therefore we can write (fTHj) as 

d /•[*! 

mq + mu 2 q± — dt'x(\t\ - t')q(±t') = -R N . (19) 

This equation can be solved for negative and positive times using the Laplace 
transformation technique. For any time dependent operator g(t), the forward 
and backward Laplace transformations are by definition 

POO 

gf{s) = I dtg(t)e- st , (20) 

and 

9 b (s) = 

/o 







g\s) = / dtg(-t)e- st , (21) 



respectively. Let x{ s ) be the Laplace transformation of x(t), then taking the 
Laplace transform of the equation (|19p. we obtain the forward and backward 
Laplace transformation of q(t), i.e., qf' b (s), as 



- m^o + ms 2 + s z x{ s ) muJ o + ms + s X\ s ) 

TsRtf(s)±R N (0) 



mu 2 + ms 2 + s 2 x(s) 



(22) 



where the upper (lower) sign corresponds to forward (backward) Laplace 
transformation respectively. 

2.2.1 A special case 

Let us as a special case, take the coupling function f(co) as 

where j3, is a positive constant. According to (|16|). this choice corresponds 
to a step function for the susceptibility 

*<H^<°; (24) 
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and accordingly, the equation (""j*?""*) can be written as 



mq + muJ q ± j3q — £(t), 



4tt 2 ,/ ^ 



e^l 



(25) 



A complete solution of this equation for t > is as follows 

q(t) = e~~%™{ — - sin u\t + g(0) coswit H — — q(0) sinc^t 

mu>i 2moJi 

R(0) /3M(Q) s M(0) . , v 

1 sm wit - v 1 sin wit - M (0) cos u) X t — sin utf] + M(t), 



m^Ml^^f^+CX*, (26) 



iii V 47T 2 7 A /u^ u; 2 - ul - &uz 



where u)\ = \J uj^ — This solution contains two parts, the first part is an 
exponentially decreasing function of time which is in fact the solution of the 
homogeneous part of the equation (|25j). the second part is a non-decaying 
term M(t), which is the response to the noise force £(£). The part M(t), 
does not have a classical correspondence and is necessary for consistency of 
the quantization of the dissipative systems. From (j""jj), (fT3J). (124*)) and (|*2*o)l . 
we can obtain an asymptotic answer for the conjugate momentum p in the 
large-time limit as 

p{t) = mM(t) + R N (t) + (3M(t) - Pq(Q). (27) 

Let the state of the damped harmonic oscillator at t = be | -0(0) ) = |0) s<S>|n) 
where |0) b, is the vacuum state of the environment and \n) is an excited state 
of the oscillator Hamiltonian, i.e., H s = 2— + ima; 2 ^ 2 , then it is clear that 

{mi ^~ + ^oV(0)|V>(0)> = (n+ \)hu Q . (28) 
From and ([?7)). we find 

*im^oo[(^(0)| : img 2 + ^mu 2 q 2 : |^(0)>] 

= [<V(0)| : ^mM 2 + ^mw 2 M 2 : |^(0)>] = 0, (29) 
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where : :, is the normal ordering operator. Therefore, macroscopically, anal- 
ogous to classical mechanics, the total energy of the particle tends to zero as 
time tends to infinity. The expectation value of the environment Hamiltonian 
(Hb), in the same state (^(O))), and large-time limit, can be calculated by 
residue calculus as 

^^[(^(0)1:^(^:^(0))] 

hpu 1 f u 2 + x 2 1 

( n +n) ~, o oxo . m n dX = {n+-)flUJ . 



nm T Jo (cj 2 - x 2 ) 2 + ^x 2 2' 

(30) 



So the total energy of the particle is absorbed by the environment and which 
is a good signal of the consistency of the approach. 

2.3 Transition probabilities 

The Hamiltonian (|%J). for a damped harmonic oscillator can be written as 

H = H + H', 

H = (a) a + -)hu + H B , 

H' = ——R+^— (31) 
m 2m' 

where a and a) are annihilation and creation operators of the Harmonic 
oscillator respectively. In interaction picture, we have 

ai (t) = e^ Hot a(0)e~^ Hot = a(0) e -^ ot , 

6 & (f) = e^\(0)e-^ = b^e**, (32) 

the terms ^p and are of the first and second order of damping respec- 
tively, therefore, for a sufficiently weak damping, ^~ is small in comparison 
with —p. Furthermore, in the first order perturbation, J^, has not any role 
in those transition probabilities where initial and final states of the harmonic 
oscillator are different from each other, hence we can neglect the term ^- in 
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H' . Substituting aj and b%j from (jB2} in -fp and using the rotating wave 
approximation [HI], one can obtain H\ in interaction picture as 



(33) 



In interaction picture, the density matrix p(t), can be obtained from the 
initial density matrix p(to) as 



Pl {t) = U I {t,t )p I {t )U}{t,t ), (34) 

where £//, is the time evolution operator, which in the first order perturbation 
is 

rt 



Uj(t,t Q = 0) « l-~ / dhH'jih), 

n Jo 



- fW4(0)e^]^. (35) 

2 

Using the density operator p/, we can calculate some transition probabilities, 
for example let pi(0) = \n)(n\ <8> |0)b b(0|, where |0)s, is the vacuum state 
of the reservoir and \n), an exited state of the harmonic oscillator, then sub- 
stituting Uj(t, 0) from ()35|) in (|34|) and tracing out the environment degrees 
of freedom, we find 



p sI {t): = Tr B ( Pl (t)), 



2 (ofr-aig) 



n^; n /-+ 00 sin z v p "V 

(36) 

where we have used the formula Ttb[\]-z)b b(1#|] — 5{k — k'). In large-time 



■ 2 ("p~"0). 
Sin — ^ 1 



limit, we can write — u>~-J ~ 2Trt6(u)f— cu ), which leads to the following 

( 2 ) 

relation for the density matrix 

(A I \/ 1 , 47T 2 ^ 3 nt|/(q; )| 2 
Psiyt) = \n){n\ H — r-^ |n-l)(n-l|. (37) 
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From p s i we can find The probability of transition |n) — > |n — 1) as 
r n ^ n _ 1 = Tr[(|n-l)(n-l|)p(t)] 

= Tr s [(\n - 1)„ w (n - l|)p s j(t)j = ^3 , (38) 

where Tr s , denotes taking trace over harmonic oscillator eigenstates. For the 
special choice (J23J), the transition probability (|3*8*Jl becomes 



ndt , . 

m 

In this case there is not any transition from the state \n) to state |n + 1). 
Now consider the case where the environment is an excited state in t = 0, 
for example, take the density matrix as 

= \ n )( n \ ® IVi'-V^s B<lft,..-lft|, (40) 

where denotes a state of the environment which contains j 

quanta with the corresponding momenta pi, ...p}, then using the relations 

Tr B [b^...l pl ) B B (l fl ,..A Pl \bl^j2 5 fi-Pi)S(k'-pi), (41) 
we find in the large-time limit, 

(-„ , IV, 3 - 2 (^Pi-^o) y. 

M') = |n)(n| + + l)<n + 1|E l/(^)l 2 ( ^-j )2 

+ ^L|n-l)(n-l| rf^zSz^, (42) 

2mh \ /\ ly^ im ^f-uoy ' v ; 

from which the transition probabilities \n) — > |n — 1) and \n) —>■ \n + 1), can 
be obtained as 

= Tr s [\n -l){n- l\p sI (t)} = -—^\f(u )\\ 
r n _ n+1 = Tr s [in + l)(n+lip s/ (t)] = ( " + J^ I/MI 2 E^-^)- 



(43) 
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Specially, for the choice (J2HJ), we have 

nf3t 



n— *n— 1 



in 



r — +1= X>("a-"o). (44) 

Another important case is when the environment has a canonical thermal 
distribution, i.e., 



6 A'T 

p 7 (0) = |n)(n| <8>p'£, = -=557, (45) 

in this case, by applying the following relations 



Trsfaplbgj = Tr B lbtp T B bt] = 0, 
e^T — l 

IV-P^l = (46) 
eA-T — 1 

we find the density operator p s j(t) in interaction picture as 
p s/ (t) := Tr B [p/(t)] 

= \n)(n\ + . ° |n+ l)(w + 1| / d 3 fc 'v 



nwo, lW /•+ 00 \j[uJr\?eT^ sin 2 ^ " o) t 
+ ^-\ n - 1 )( n ~ 1 \ d 3 k^\Jf -— 2 — , (47) 



2mh' " '7-oo e ^_x (53=^)2 

which accordingly lead to the following transition probabilities in large-time 
limit 

f T n 1\/ ii An2uJ o nt l/N| 2 e^ 

r n ^ n _i = Tr s [\n - l)(n - l|p fl /(t)J - 



r n ^„ + i = Tr s [|n + l)(n+ 1| 



mhc 3 e"^r — 1 
47r^ 3 (n + l)t|/(a; )| 2 



mrt-C 5 g ITT — 1 

(48) 
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and for the special case are reduced to 



r 



n(3tei<T 




1) 



r 



(n + l)0t 



(49) 




1) 



Therefore, in low temperatures regime, the energy flows from the oscillator 
to the environment by the rate r n _» n _i i— > 22, and no energy flows from the 
environment to the oscillator in this case. 

3 Three-dimensional quantum dissipative sys- 
tems 

In this section, we study a general three-dimensional quantum dissipative 
system. This time for generality, we model the environment of the system by 
two independent quantum fields, namely B and B quantum fields. Quantum 
field B, interacts with the momentum of the main system through a minimal 
coupling term and the quantum field B, interacts with the position operator 
of the main system similar to a dipole interaction term. For this purpose we 
take the Hamiltonian of the environment as 



where according to the explanations of the previous section, we choose = 



H E = H B (t) + Hz(t) 




(50) 



c\k\. The annihilation and creation operators 
commutation relations 



&£, bp and at, satisfy the 





(51) 



and the rest of commutation relations are zero 
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Let the damped system be a particle with mass m under an external po- 
tential v(x), we take the total Hamiltonian, i.e., system plus the environment, 

as 

H= iP-R) +v (£)-R.£ + h b + H § , (52) 

where x and p, are position and canonical conjugate momentum operators of 
the particle respectively, and satisfy the canonical commutation rules 

[xi,pj] = ihSij. (53) 

Operators R and R, play the basic role in interaction between the system 
and the environment and are defined by 

-, r+oo 

R(t) = / d'kifi^it) + F^bim, 

J-oo K 

R(t) = / d s k[g(u^(t) + g^dlm, (54) 

J — oo 

where k — 4r. For future purposes the fields R and R, are taken to have 
|fc| 

the longitudinal polarization, et us call the functions f{oj£) and g{u^j:), the 
coupling functions between the environment and the system. The Heisenberg 
equations for the position x(t) and the momentum p of the particle are 

x = =[H,x\ = , 

n m 

p= ~[H,fl = -Vv + R. (55) 

By eliminating p in relations ()55|h we come to the following equation for the 
damped system 

mS = -Vf — R + R. (56) 

Using (|51|). one can easily find the Heisenberg equations for the operators 
and d% as 

z z " 

d * = n^' = + n 9 *^k) k ■ f • ( 57 ) 
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These equations can be solved formally as 



(58) 



Substituting these solutions into (f51J). one obtains what we have called the 
constitutive equations of the environment, i.e., 

R(t) = R N (t) + / dt' X (\t\ - t')x(±t'), 
Jo 

* * r\t\ 

R(t) = R N (t) + / dt'x(\t\ - t')x(±t'), (59) 
Jo 

where the upper(lower) sign, corresponds to t > 0(t < 0), respectively. The 
functions 

167T f°° 

X ^ = 3h^Jo rfu;w2 |/M| 2 sin^ t>0, 

X (t) = t<0, (60) 



and 



167T f°° 

X(t) = J o duuj 2 \g(uj)\ 2 sin ut t > 0, 

m = o t<o, (6i) 

are called the susceptibilities of the environment. The operators Rn and Rn, 
are the noise operators 

-» f+oo 

R N (t) = / d 3 M/H)%(0)e-^ + /* K)4 (0)e^*]fc, (62) 

J-oo re 
2 f+oo 

i2jv(*) = / d 3 fe^(w £ )d jE (0)e- w s' + g*(uM(0)e^]k. (63) 

J-oo " K 

If we are given the macroscopic susceptibilities x(£) and x(t), which are zero 
for t < 0, then we can invert the equations ()60|) and ()61|) and obtain the 
corresponding coupling functions f(u>) and g>(o;) as 



driC r°° 

\f(uo)\ 2 = ——^ / dt X (t)sinut, u>0, 
8ty z uj z Jo 

l/MI 2 = 0, w = 0. (64) 
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I^MI 2 = 022 dtx(t)sinut, uj>0, 
8n z u z Jo 

\g{u)\ 2 = 0, u = 0. (65) 

It should be noted that the sign of the right hand side of the formulas (|64|) 
and should be positive, this is because, the right hand side of these for- 
mulas is proportional to the imaginary part of the susceptibilities x(i) and 
x(t) in the frequency domain and this imaginary part is also connected to 
energy losses of a system in the presence of an environment, so it is necessar- 
ily positive. If these integrals would not be positive for all frequencies, the 
susceptibilities x(t) and x(t) should be discarded, as they would be unphys- 

ical. Substituting R and R, from the constitutive equations (}5T?|) . into (j5EJ), 
one can obtain the equation of motion of the operator x as 

d r\M ■ r\*\ -> u ? 

mx± — dt'x(\t\-t')x(±t')- dt'x(\t\-t')x(±t') = -Vv(x)-R N +R N , 
at Jo Jo 

(66) 

where the upper (lower) sign corresponds to t > (t < 0). This is the gener- 
alized macroscopic Langevin-Schrodinger equation for a dissipative quantum 
system. Therefore, coupling with the environment in microscopic level, is 

described macroscopically by two random forces —Rn and Rn and the cor- 
responding memory functions x an d X- 

As a simple example, let us consider a damping three-dimensional har- 
monic oscillator with resonance frequency oj and mass rn, where the damping 
force is proportional to the velocity and position of the oscillator. In this case, 
we take the susceptibilities %(t) and x(i), as follows 

x(t) = 





< t < A, 
otherwise, 



where < a < 1 and (3 are some positive constants. Now one can use the 
relations (J64|) and (J65|) to obtain the related coupling functions as 

I/(W)|2 = S0' (67) 
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. ,, 9 3hc 3 amuj?: sin 2 ^ 

= SirV f ' ^ (68) 



In this case, the Langevin-Schrddinger equation (jobj) . is reduced to 

amw 2 /■!*! 



mx + mu x ± px — / at x(±t ) = —Rn + Rn, 

A 7|t|-A 

(69) 



where Rn and R are the noise operators (j62j) and (|63|) with coupling functions 
()67|) and ()68|) respectively. It is clear that in the limit A — ► 0, the coupling 

function (}68|) . tends to zero and accordingly, the noise force Rn vanishes. 
Therefore, the equation (jfinj) becomes 



mx + (1 — a)nvjj\x ± {3x = £(t), 




(70) 

For an absorptive environment, < a < 1, otherwise the environment would 
behave as an amplifier, i.e., energy transfers to the system from the envi- 
ronment, this may be an indication for making other related models, but 
in our discussion, this case is not physical. As is seen from this equation, 
the environment affect the motion of the particle by two damping forces, 
one is a friction force proportional to the velocity and the other is a force 
proportional to the position and in opposite direction of the external force 
— Vv(x) = — mijj\x. In fact, the later damping force, influence the frequency 
of the oscillator. 



3.1 Dissipative quantum two- level systems 

The rate of spontaneous decay of a two-level quantum system imbedded in an 
absorptive environment, can be calculated in this approach as follows. Here 
we assume that the interaction between the particle and its environment is so 
weak that we can use the Weiskopf-Wigner approximation for calculating 
the decay constant. Let us write the total Hamiltonian (|52j). as the following 
form 

H = H + H', 
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H = ^- + v{x) +H B + H s , 
2m 

H' = -—■ R-R-x, (71) 
m 

where we have ignored from the term J^-, because, this term does not con- 
tribute in the following calculations of the decaying rate of the system. Now 
for a considerable simplification, we restrict ourselves to a two-state model 
of a system. In fact, the Hilbert space of the main system is artificially trun- 
cated to the two states |1) and |2) with unperturbed energy eigenvalues E x , 
E2, respectively. Therefore, for the Hamiltonian of the system we can write 

H s = ^ + v{x) = hn a z , a z = |2)(2| - |1)(1|, fi = E2 ~ E \ 
2m 2 n 

(72) 

and the interaction Hamiltonian H' becomes 

H' = -i{l J d*k[ab z f(u % )x l2 ■ k + abif*(^)x 12 ■ k - C.C] 

- J d 3 k[ad^g(u^)x 12 ■ k + adtg*^)^ ■ k + C.C]., (73) 

where 

a=\l}(2\, x 12 = (l|f|2), (74) 

and we have used this fact that the energy eigenstates of the Hamiltonian 
(172^1 have a well defined parity, so that the diagonal elements of the operators 
x and p, are zero. 

To study the spontaneous decay of an initially excited two-level system, 
imbedded in an absorptive environment, we may look for the system wave 
function at time t in interaction picture 



= c(t)\2)\0) B \0) B + J d 3 kD,:(t)\l)\k) B \0) B + 
d 3 kG % (t)\l)\0) B \k) B , (75) 



where \0) B and \0) B , are vacuum states of quantum fields B and B, respec- 
tively, and the coefficient c(t), Dg(t) and G^(t), are to be specified by the 
Schrodinger equation 

lh dT^ = H i® = ^H' s e^, (76) 



19 



with initial conditions c(0) = 1, D%(0) = Gg(0) = 0. Where Hq and H' s are 
Hamiltonians Hq and H' in the Schrodinger picture respectively. Substituting 
\if}(t))i from (J7SJ), in (|75jl. and using the rotating wave approximation 
we find the differential equations 



l fcc(i) = -/ d 3 A;[-zfio^*i2 ■ kfiu^e^-^D^t) 

+ ^■kg^y^-^G^t)}, (77) 

zftl^i) = -A>£i2 ■ fcr(^)e- l(n °-^ c (t), (78) 

aftG^t) = _f 12 . fc/(c^) e -^°-^) f c (t), (79) 

One can solve this differential equations using the Laplace transformation 
technique and deduce 

c(t) = C dt'7(t-t')c(t'), 
Jo 

j(t-t') = -J d 3 k[nl\x u -kf(uj:)\ 2 + \x u -k^^ 

(80) 

Here we restrict our attention to the weak coupling regime, where the Markov 
approximation [36J applies. That is to say, we may replace c(t') in the inte- 
grand of (|80|) with c(t) and approximate the time integral f dt'j(t — t'), in 
large-time limit, as 



dtle i[n -^)(t-t') w [lP + n5 ^ Q _ ^ (gl) 

' '0 ~~ 

where P, denotes th principal Cauchy value. After some algebraic calcula- 
tions, we find 

c(t) = -(P + iA)c(t), (82) 
where (3 is the decay constant 



_ 47r 2 f^i2| 2 |/(f2o)| 2 47r 2 |f 12 | 2 |^(fi ) 

P ~ o+2 Q + 



3ft 2 c 3 3h 2 c 3 



(83) 



and A is the shift frequency 



3ft 2 c 3 7o Vln-u V 7 
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4 Dissipative scalar field theory 



In this section, we discuss the dissipative scalar field theory. We take for 
convenience a 1+1-dimensional scalar field but the generalization to a general 
scalar field theory is straightforward. For this purpose, assume ip(x,t) be a 
scalar field operator defined on a closed compact interval [0, L], satisfying the 
boundary conditions ^(O) = ifi{L) = 0. The field ip(x,t), can be expanded in 
terms of the orthogonal wave function sin as 



il>(x,t) = E \lj^K(t) + 4(t)]sin^, (85) 



n=l 



where uj n = yf^, and X, fi, are some constants. The operators a n and 
are annihilation and creation operators of the scalar field ip(x,t), respectively 
and satisfy the following equal-time commutation relations 

\a n (t),a ] m (t)] = S nm . (86) 

The conjugate canonical momentum density of the field ip(x,t), can be also 
expanded in the same basis as 



CO 



7T^(x,t) = I 2^ y — —[<(*) ~ a n{t)\ S1I1_ ^- ( 87 ) 

From (|HU|) . we deduce that ip, tt^,, satisfy the equal time commutation relation 

[ifj(x,t),ii^(x' ,t)} = ih5(x — x'). (88) 
The Hamiltonian of the scalar field ifj(x,t), in normal ordered form is 

. r id i 

n=l 



^ = / o : + := E ^n^KW, ( 89 ) 



where i/^, denotes the derivative with respect to x. Now similar to what 
we did in previous sections, here we model the environment by a quantum 
field B containing an infinite number, but numerable, of Klein-Gordon fields. 
Therefore, the environment Hamiltonian, can be written as 

H B (t) = E / Mtuoibyfibjfi, w s = c\k\, (90) 

n=l J -°° 
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where b n ^ and p are annihilation and creation operators of the environment 
respectively and satisfy the following bosonic commutation relations (bosonic 
bath) 

[te(*),&U(*)] = <w(£-#), (91) 

In a systematic way, and also according to what we did in the previous 
sections, we write the total Hamiltonian, i.e., the scalar field ip(x,t) plus the 
environment, as 

= f L dx M^)-R(x,t) ? + l + ^ 



Jo ~~ 2A ' 2^ rx 

wherein the operator R(x,t), has the basic role in interaction between the 
scalar field and its environment and is defined by 

/2 00 r+oo , tittt 

iT,]^ rf 3 M/K^)^W + fk^)^W]sm T . (93) 

Let us call the function f(uk,x), the coupling function between the scalar 
field and its environment. The coupling function is independent (dependent) 
on x for a homogeneous(inhomogeneous) environment respectively. 

One can easily show that the Heisenberg equations for ip(x, t) and 7Tf(x, t) 

are 

• , . 1 rTT ip(x,t) — R(x,t) 
if;{x,t) = ^[H, ifi] = ^ , 

I 

~Kt/i(x,t) = —[HjTTip] — —fJ>lf>xxi (94) 

we can eliminate between the relations (|94|) and finally find the following 
equation for the field if} 

Xip - iiip xx = -R. (95) 

Now similar to the previous sections, if we write the Heisenberg equation 
for the annihilation operator 6 n g and solve it formally then substitute this 
solution into (J93j) . we can obtain a constitutive equation for the environment 

as 

R(x, t) = R N (x, t) + / dt'x(x, \t\ - t')ip(x, it 7 ), (96) 



where the upper(lower) sign, corresponds to t > 0(t < 0) respectively and 



57T 



x( x ,t) = ^7 / duu \f(u),x)\ sina;t, t > 0, 

he 6 Jo 

X (x,t) = 0, t<0, (97) 
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is the susceptibility of the environment. By inverting the relation (J97|) . we 
can obtain the coupling function f(uj,x), in terms of the susceptibility, as 

\f(u,x)\ 2 = 4 ^ 2 ^ 2 J dtx(x,t)smut, u>0, 
\f(u,x)\ 2 = 0, u = 0. (98) 

The operator 

jY. loo d M/(^^)^(0)e-^+rK,x)^( )e^]sin — 

(99) 

is a noise operator and necessary for a consistent quantization of a dissipative 
theory. Finally, substituting (f9T)J) in (jHSJ), leads to the following wave equation 
for the field operator -0 

Xip - iiip xx ± — / dt'x(x, \t\ - t)ip(x, ±t) = -R N (x,t). (100) 
dt Jo 

Taking the Laplace transformation of this equation, we find that the forward 
and backward Laplace transforms of ip(x,t), satisfy the following equation 

[As 2 + s 2 x(x, s)]t f '\x, s) - ^ b x (x, s) = Jf> b (x, s) (101) 

where 



jf> b ( X} s ) = [As + sx(x, s)]i)(x, 0) =f Xtp(x, 0) =f Rtf{x, s) ± R N (x, 0), 

the upper(lower) sign, corresponds to the forward (backward) Laplace 
transformations respectively and x{ x i s )-> ls the Laplace transform of the sus- 
ceptibility. One can solve equation ()101|) in terms of the Green function 

^x, s ) = [ L dx'G(x, x', s)J f > b (x', s), (102) 
— Jo 

where G(x,x',s) satisfies 

(As 2 + s 2 x(s, x))G(x, x' , s) — fj,G xx (x, x', s) = S(x — x'), (103) 

with the boundary conditions G(0,x',s) = G(L,x',s) = 0. For a homo- 
geneous environment, the coupling function and the susceptibility x-, are 
position independent and in this case, the Green function is 

2 ^ 1 . nirx . nnx' 

G{x, x , s = - ^ sm — — sin — — . 104 

L As 2 + s 2 x(s) + \uji L L 
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4.1 A special case 

For the special choice of coupling function (J23J), equation (jlOUj) . becomes 
A-0 - fiipxx ±f3ip = i{x, t), 



*) = «V^37 E / -=(U0)e-^ - ^(O)e^] sin 



L 
(105) 



By defining 



Xn = V^T^" + Pn = i\J T ^Y L {al-a n ), (106) 

and using (JTU3J), we obtain 

= i\h4r 2 / -7=y0)e^ - &i,(0)e*n (107) 



a complete solution of this equation for t > 0, is 
x n (t) = e~^{x n (0) cosQ n t + 77777-^.(0) sinf] w t - sin ^ 
- M n (0) cosf^ + Xw(0) ~ Mw(0) sinf^} + M n (t), 



r+00 

M n {t) = i I d 3 k A 



^ r 6^(0) ^-(0) rtorf) 



N ^^^^ ^ - "| - if^S ^ - W| + 

(108 



where f2 n = ycj 2 — This solution consists of two parts, the first part 
is an exponentially decreasing function in time which is the solution of the 
homogeneous part of the equation ([107|h i.e., when ( n (t) = 0, the second part 
is M n (t), which is an oscillatory exponential function and is the response to 
the noise force Cn(t)- The answer M n (t), does not have a classical counterpart 
and is necessary for a consistent treatment of dissipative quantum systems. 
In fact, without M n (t), the commutation relations [x n ,p n ], tend to zero in 
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large-time limit and accordingly, the uncertainty relations would be violated. 
From (J93J), and (jl(J8|) . we obtain the asymptotic answer 

Pn (t) = XM n (t) + r N (t) + (3M n (t) - Px n (0), 
r N (t) = E /_ m 7=[^(0)e-^ + ^(0)e*n (109) 



V fe 

for Pn(t), when t — > +00. If the state of the system in t — 0, is taken to be 
|^(0)) = |0) B ® |l mi , l mr ) s , where |l mi , l mr ) s is an excited state of the 
Hamiltonian H s , then it is clear that 

00 r 

(V(0) E na;„at (0)a n (0)|^(0)) = E^> ( 110 ) 

n=l i=l 

on the other hand, from (jlOSJ) . we find 



Um^imi ■ f L dx[\\i, 2 + : |^(0))] 

Jo 2 I 

00 1 1 

MO) I : J2(n XM n + ~M n Ml) : |^(0)) = 0. 

n=l Z Z 



111) 



So in the large-time limit, the expectation value of total energy of the system 
in the initial state |^>(0)) = |0)# <8> |l mi , l mr )s, tends to zero, as expected. 
Also, using (|109|) . it is clear that 

Um^Kmi : E ff + = 1^(0))] = fr<4(«)> = ^2 E 7— ■ 

(112) 

Now, if one obtains the Heisenberg equation for b n %(t) and solve it formally 
in terms of x n (t), then using (jl(J8|) . he (she) will find 



limt^wm ■■ / AEH^W^W : l^(o))] 

n=l 



ftp v- r 

7tA i=1 Jo 



r 



dx, 



=1 



E^mi- (113) 
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Therefore, the energy of the system is completely transferred to the environ- 
ment, as expected macroscopically. 



4.2 Transition probabilities 

As in the second section, we can calculate some of transition probabilities 
for the scalar field ip by tracing out the environment degrees of freedom. For 
example, Let the total density operator in the interaction picture at time 
t = 0, be 

Mo) = I I'm) s s(j"m |®|0)bb<0|, (114) 

where |0) B , is the vacuum state of the environment and \r m ) s = ^^-\0) s is 
an excited state (Fock state) of the field tp, in large-time limit, we obtain the 
following probability transitions 

r |r m > s -|r m -i> s = Tr[\r m -l) s s {r m - l\p(t)}, 

= Tr s [\r m -1) 8 a (r m -l\p aI (t)], 

where Tr s , means taking trace over the scalar field eigenstates. 

Another important case is when the environment has a canonical thermal 
distribution 

-Hg 

P/(0) = \r m ) s s(r m \ ® Pb, Pb = (H6) 



TR B {e 

then in large-time limit, we obtain the following transition probabilities 

4n 2 ujlrt\f{u m )\ 2 ei» 

\r m )s^\r m -l)s Xhc3 e f^ _ 1 ' 

Xhc 6 e kt — 1 

{T + l)^ m t\f{u m )\ 2 
t\r m ) s -+\r m+ l) s ^3 e ^_l 

So in low temperatures, the energy flows from the oscillator to the it's en- 
vironment by the rate ^ and no energy flows from the environment to the 
scalar field. 
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5 Dissipative vector field theory 



Quantization of a quantum vector field Y in a three-dimensional inhomoge- 
neous absorptive environment can be investigated by modeling the environ- 

— * — 

ment of Y by two independent quantum fields, namely B and B quantum 
fields. The susceptibility of the environment and the quantum noise fields are 
identified in terms of ladder operators of the environment and parameters of 
this model as in previous sections. We assume that the vector field Y can be 
propagated in infinite space with a suitable boundary condition at infinity, 
that is, the field Y tends to zero at infinity. In this case both the fields B 
and B, contain a continuum of Klein-Gordon fields. It is also remarkable to 
note that if the volume in which the field Y can be propagated is a finite 
volume, for example a cubic cavity, then the fields B and B, will contain a 
numerable set of Klein-Gordon fields, as in the previous section for the case 
of a scalar field. 

From the interaction point of view, the field B interacts with the con- 
jugate canonical momentum density of the main vector field Y through a 
minimal coupling term and quantum field B interacts with the field Y sim- 
ilar to a dipole interaction term. In this scheme, we take the environment 
Hamiltonian as 



He = Hb + H 



B ' 

3 



H B = S tf d 3 k f d 3 qhu,;bl(k,q,t)b u (k,q,t), 

u=l J J 

3 

H B = J2 I d3 k I ^qhu^dKk, q, t)d v (k, q, t). (118) 

where the annihilation and creation operators b„(k, q, t), b\,{k, q, t), d\,{k, q, t) 
and dl(k,q,t) of the environment, satisfy the commutation relations 

[K(k, q, t),bl(9, q', t)\ = 5 uu ,5(k - k')8(q- q'), 
[d u (k, q, t),dl(k>, q', t)} = 5 vv ,5{k - k')S(q- q<). 

(119) 

Let us assume that the Hamiltonian of the main system is 

H Y = J d^ + l -p{r)ul{r)Y\ (120) 
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where 7?y, is the conjugate canonical momentum density of Y. An example 
of this kind of Hamiltonian, the vector field Y, can be an electric polarization 
density in an absorptive dielectric medium with related eigenfrequency u)o(r) 
and density p{r) |37j. 

The conjugate fields Y and satisfy the commutation relation 

\Y i (r,t),ti Yj (?,t)]=ih5 ij 6(f-P). (121) 

If the damping forces together with the restoring force —p(f)ujQ(r)Y , are 
exerted on the elements of the medium, then the equation of motion of Y 
should be a Langevin-Schrodinger equation (J66|) . wherein, the external force 
— Vi>, is replaced with — p{f)ul{r)Y . For this purpose, we take the total 
Hamiltonian, i.e., the main system (vector field Y), plus the environment as 

H = j d 3 ^^ 2 + \l*f)<4lffi* — R-Y + H B + Hg. (122) 

The operators R(r,t) and R(f,t), play the basic roles in the interaction 
between the environment and the system and are defined by 

R(r,t) = £ f & I -P=\}{^K$M<!** 
v=i J J J (2tt) 3 

(123) 

R(r,t) = ±[d 3 k[ ^J=[g(^r)d^q,t)e^ 
u=i J J y(27r) 3 

(124) 

where 



u= 1,2, 



u 3 (q) = q= —, 
m 

q-e vq *= 0, 



^ = 1,2, 



(125) 
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are three orthogonal unit vectors for any q. 

The function f(u^,r) and f{u!^, r), are the coupling functions, which are 
position dependent (independent) for an inhomogeneous(homogeneous) envi- 
ronment respectively. The equations of motion for the fields Y and fry, can 
be obtained from the Heisenberg equations 

dt h [ ' 1 P ' 

d -^ = \[H^ Y } = -pul{r)Y + R. (126) 
By eliminating fry between these equations we obtain the equation of motion 

— * 

of the vector field Y as, 

4 + /W » (rT = T +A (127) 

Similar to what we did in the previous sections one can obtain, the constitu- 
tive equation of the environment as 



R(f, t) = R N (f, t) + / dt' X (f, \t\ - t')Y(f, ±t'), 
Jo 

R(f, t) = R N (f, t) + /"'*' dt'x(r, \t\ - t')Y(f, ±t'), (128) 
Jo 



where the upper (lower) sign corresponds to t > 0(t < 0), respectively. The 
following relations 

X (r,t) = t<0 (129) 

gyj- r-oo 

x(^0 = ^y o duu 2 \g(u,r)\ 2 smujt, t > 0, 

X(r,t) = 0, t<0, (130) 

give the susceptibilities of the environment in terms of the coupling functions. 

— * — 

Operators Rjkj and R N , are noise fields 



3 ,g _» 

R N (r,t) = j d*k j -j=L=[f{u % ,^K(k ,q,0)e 
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+ r(co % , r)bt(k, q, 0)e^-^%(g) 



(131) 



R N (f,t) = Yl d 3 k / -^[g^^d^q^e 




+ g*{u %1 r)dl(k, q, ti)e^-^]u v {q)- 



(132) 



Finally, substituting ()128f) in ()127|) . we find a generalized Langevin-Shrodinger 



This equation can be solved using the Laplace transformation technique 
for negative and positive times similar to the previous sections. 

6 Concluding remarks 

By considering a dissipative quantum system as an open system which in- 
teracts with it's environment, a modeling of the environment by massless 
Klein-Gordon fields, is achieved. Using the above idea, i.e., the field version 
of an environment, an arbitrary dissipative quantum system is quantized sys- 
tematically and consistently. Some coupling functions are introduced which 
physically describe the environment under consideration and also have a ba- 
sic role in interaction between the system and it's environment. Inspired by 
the ideas of electrodynamics in a media, some susceptibility functions are at- 
tributed to the environment and formulas connecting these susceptibilities to 
the coupling function, are found. The quantum Langevin-Schrodinger equa- 
tion is obtained directly from the Heisenberg equations. The explicit form 
of the noise terms are obtained. Some transition probabilities indicating the 
way energy flows from the system to it's environment, are calculated and the 
energy conservation, as a signal of consistency, is explicitly worked out. The 
whole formalism is generalized to the case of a dissipative scalar and vector 
field theory straightforwardly. 



equation for the damped vector field Y 



pY + pu 2 (r)Y ± j t jf dt' X (f, \t\ - t')Y(r, ±t') - 
[ t] dt'x(f, \t\ - t')Y(f, ±t') +R N - R N £(r, t) = 



(133) 
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